Landscape genomics reveals regions associated with adaptive phenotypic and genetic variation in Ethiopian indigenous chickens.

Climate change is a threat to sustainable livestock production and livelihoods in the tropics. It has adverse impacts on feed and water availability, disease prevalence, production, environmental temperature, and biodiversity. Unravelling the drivers of local adaptation and understanding the underlying genetic variation in random mating indigenous livestock populations informs the design of genetic improvement programmes that aim to increase productivity and resilience. In the present study, we combined environmental, genomic, and phenotypic information of Ethiopian indigenous chickens to investigate their environmental adaptability. Through a hybrid sampling strategy, we captured wide biological and ecological variabilities across the country. Our environmental dataset comprised mean values of 34 climatic, vegetation and soil variables collected over a thirty-year period for 260 geolocations. Our biological dataset included whole genome sequences and quantitative measurements (on eight traits) from 513 individuals, representing 26 chicken populations spread along 4 elevational gradients (6-7 populations per gradient). We performed signatures of selection analyses ([Formula: see text] and XP-EHH) to detect footprints of natural selection, and redundancy analyses (RDA) to determine genotype-environment and genotype-phenotype-associations. RDA identified 1909 outlier SNPs linked with six environmental predictors, which have the highest contributions as ecological drivers of adaptive phenotypic variation. The same method detected 2430 outlier SNPs that are associated with five traits. A large overlap has been observed between signatures of selection identified by[Formula: see text]and XP-EHH showing that both methods target similar selective sweep regions. Average genetic differences measured by [Formula: see text] are low between gradients, but XP-EHH signals are the strongest between agroecologies. Genes in the calcium signalling pathway, those associated with the hypoxia-inducible factor (HIF) transcription factors, and sports performance (GALNTL6) are under selection in high-altitude populations. Our study underscores the relevance of landscape genomics as a powerful interdisciplinary approach to dissect adaptive phenotypic and genetic variation in random mating indigenous livestock populations.

The origin of genetic diversity of indigenous cockfighting chickens of Pakistan by analyzing the mtDNA.

In Pakistan, the origin of the indigenous cockfighting chicken (ICC) or gamecock population is unknown. However, it is speculated that this might have been associated with domestication, an event linked to recreational, entertainment (cockfighting), religious or ornamental activities. This study aims to understand the origin and genetic diversity of the ICC population in Pakistan. A total of 185 ICC population and 10 captive Indian red junglefowl (Gallus gallus murghi) were analyzed for genetic diversity indices and phylogenetic reconstruction using a 397 bp of mtDNA D-loop region. It is reported that a total of 43 haplotypes from 38 polymorphic nucleotide sites. The haplotype and nucleotide diversity are also estimated in the range of 0.643-0.909, and 0.00585-0.01575, respectively. The total genetic diversity within the population was 91.52%. Four mitochondrial haplogroups A, B, C and D were identified by median-joining network analysis, two of them have high percentages, haplogroup D (81.6%) and A (15.1%). Phylogenetic analysis showed that the ICC population of Pakistan and Gallus gallus murghi shared haplogroup D. The results of this study showed that sub-haplogroup D17a05, has significantly high haplotype diversity and percentage as compared to previously published studies, this indicated that Pakistan might be one of the centres of domestication for chicken, as it is considered that Southeast Asia is the centre of domestication. Frequencies of Haplogroup A also indicate South-North indices. This research work showed that the indigenous cockfighting chicken population of Pakistan is genetically introgressed from Gallus gallus murghi, and significant variations could be attributed to the underlying differences in the geographics, selection pressures, introgression, and regional practices; and multiple origins of cockfighting chickens' populations around the world which reflected the past trading routes between human communities and civilizations.

Exploring genetic diversity and variation of Ovar-DRB1 gene in Sudan Desert Sheep using targeted next-generation sequencing.

INTRODUCTION: The Ovar-DRB1 gene, a crucial element of the Major Histocompatibility Complex (MHC) Class II region, initiates adaptive immunity by presenting antigens to T-cells. Genetic diversity in sheep, particularly in MHC Class II genes like Ovar-DRB1, directly influences the specturm of presented antigens impacting immune responses and disease susceptability. Understanding the allelic diversity of Ovar-DRB1 gene in Sudan Desert Sheep (SDS) is essential for uncovering the genetic basis of immune responses and disease resistance, given the the breeds significance in Sudan's unique environment. METHODS: Utilizing Targeted Next-Generation Sequencing (NGS) we explore allelic diversity in Ovar-DRB1 gene within SDS. Successfully ampliying and and sequencing the second exon of this gene in 288 SDS samples representing six breeds provided a comprehensive allelic profile, enabling a detalied examination of the gene's genetic makeup. RESULTS: We identifed forty-six alleles, including four previously unreported, enrichness the genetic diversity of SDS breeds. These alleles exhibiting non-uniform distribution, varying frequencies across breeds, indicating a breed-specific genetic landscape. Certain alleles, known and novel, show higher frequencies in specific populations, suggesting potential associations with adaptive immune responses. Identifying these alleles sets the stage for investigating their functional roles and implications for disease resistance. Genetic differentiation among SDS breeds, as indicated by FST values and clustering analyses, highlights a unique genetic makeup shaped by geographic and historical factors. These differentiation patterns among SDS breeds have broader implications for breed conservation and targeted breeding to enhance disease resistance in specific populations. CONCLUSION: This study unveils Ovar-DRB1 gene allelic diversity in SDS breeds through targeted NGS and genetic analyses, revealing new alleles that underscore the breeds' unique genetic profile. Insights into the genetic factors governing immune responses and disease resistance emerge, promising for optimization of breeding strategies for enhanced livestock health in Sudan's unique environment.

Genomic analysis of Nigerian indigenous chickens reveals their genetic diversity and adaptation to heat-stress.

Indigenous poultry breeds from Africa can survive in harsh tropical environments (such as long arid seasons, excessive rain and humidity, and extreme heat) and are resilient to disease challenges, but they are not productive compared to their commercial counterparts. Their adaptive characteristics are in response to natural selection or to artificial selection for production traits that have left selection signatures in the genome. Identifying these signatures of positive selection can provide insight into the genetic bases of tropical adaptations observed in indigenous poultry and thereby help to develop robust and high-performing breeds for extreme tropical climates. Here, we present the first large-scale whole-genome sequencing analysis of Nigerian indigenous chickens from different agro-climatic conditions, investigating their genetic diversity and adaptation to tropical hot climates (extreme arid and extreme humid conditions). The study shows a large extant genetic diversity but low level of population differentiation. Using different selection signature analyses, several candidate genes for adaptation were detected, especially in relation to thermotolerance and immune response (e.g., cytochrome P450 2B4-like, TSHR, HSF1, CDC37, SFTPB, HIF3A, SLC44A2, and ILF3 genes). These results have important implications for conserving valuable genetic resources and breeding improvement of chickens for thermotolerance.

Whole-genome resource sequences of 57 indigenous Ethiopian goats.

Domestic goats are distributed worldwide, with approximately 35% of the one billion world goat population occurring in Africa. Ethiopia has 52.5 million goats, ~99.9% of which are considered indigenous landraces deriving from animals introduced to the Horn of Africa in the distant past by nomadic herders. They have continued to be managed by smallholder farmers and semi-mobile pastoralists throughout the region. We report here 57 goat genomes from 12 Ethiopian goat populations sampled from different agro-climates. The data were generated through sequencing DNA samples on the Illumina NovaSeq 6000 platform at a mean depth of 9.71x and 150 bp pair-end reads. In total, ~2 terabytes of raw data were generated, and 99.8% of the clean reads mapped successfully against the goat reference genome assembly at a coverage of 99.6%. About 24.76 million SNPs were generated. These SNPs can be used to study the population structure and genome dynamics of goats at the country, regional, and global levels to shed light on the species' evolutionary trajectory.

Mapping restricted introgression across the genomes of admixed indigenous African cattle breeds.

BACKGROUND: The genomes of indigenous African cattle are composed of components with Middle Eastern (taurine) and South Asian (indicine) origins, providing a valuable model to study hybridization and to identify genetic barriers to gene flow. In this study, we analysed indigenous African cattle breeds as models of hybrid zones, considering taurine and indicine samples as ancestors. In a genomic cline analysis of whole-genome sequence data, we considered over 8 million variants from 144 animals, which allows for fine-mapping of potential genomic incompatibilities at high resolution across the genome. RESULTS: We identified several thousand variants that had significantly steep clines ('SCV') across the whole genome, indicating restricted introgression. Some of the SCV were clustered into extended regions, with the longest on chromosome 7, spanning 725 kb and including 27 genes. We found that variants with a high phenotypic impact (e.g. indels, intra-genic and missense variants) likely represent greater genetic barriers to gene flow. Furthermore, our findings provide evidence that a large proportion of breed differentiation in African cattle could be linked to genomic incompatibilities and reproductive isolation. Functional evaluation of genes with SCV suggest that mitonuclear incompatibilities and genes associated with fitness (e.g. resistance to paratuberculosis) could account for restricted gene flow in indigenous African cattle. CONCLUSIONS: To our knowledge, this is the first time genomic cline analysis has been applied to identify restricted introgression in the genomes of indigenous African cattle and the results provide extended insights into mechanisms (e.g. genomic incompatibilities) contributing to hybrid differentiation. These results have important implications for our understanding of genetic incompatibilities and reproductive isolation and provide important insights into the impact of cross-breeding cattle with the aim of producing offspring that are both hardy and productive.

Inference of Admixture Origins in Indigenous African Cattle.

Present-day African cattle retain a unique genetic profile composed of a mixture of the Bos taurus and Bos indicus populations introduced into the continent at different time periods. However, details of the admixture history and the exact origins of the source populations remain obscure. Here, we infer the source of admixture in the earliest domestic cattle in Africa, African taurine. We detect a significant contribution (up to ∼20%) from a basal taurine lineage, which might represent the now-extinct African aurochs. In addition, we show that the indicine ancestry of African cattle, although most closely related to so-far sampled North Indian indicine breeds, has a small amount of additional genetic affinity to Southeast Asian indicine breeds. Our findings support the hypothesis of aurochs introgression into African taurine and generate a novel hypothesis that the origin of indicine ancestry in Africa might be different indicine populations than the ones found in North India today.

Candidate signatures of positive selection for environmental adaptation in indigenous African cattle: A review.

Environmental adaptation traits of indigenous African cattle are increasingly being investigated to respond to the need for sustainable livestock production in the context of unpredictable climatic changes. Several studies have highlighted genomic regions under positive selection probably associated with adaptation to environmental challenges (e.g. heat stress, trypanosomiasis, tick and tick-borne diseases). However, little attention has focused on pinpointing the candidate causative variant(s) controlling the traits. This review compiled information from 22 studies on signatures of positive selection in indigenous African cattle breeds to identify regions under positive selection. We highlight some key candidate genome regions and genes of relevance to the challenges of living in extreme environments (high temperature, high altitude, high infectious disease prevalence). They include candidate genes involved in biological pathways relating to innate and adaptive immunity (e.g. BoLAs, SPAG11, IL1RL2 and GFI1B), heat stress (e.g. HSPs, SOD1 and PRLH) and hypoxia responses (e.g. BDNF and INPP4A). Notably, the highest numbers of candidate regions are found on BTA3, BTA5 and BTA7. They overlap with genes playing roles in several biological functions and pathways. These include but are not limited to growth and feed intake, cell stability, protein stability and sweat gland development. This review may further guide targeted genome studies aiming to assess the importance of candidate causative mutations, within regulatory and protein-coding genome regions, to further understand the biological mechanisms underlying African cattle's unique adaption.

Genome-wide diversity and admixture of five indigenous cattle populations from the Tigray region of northern Ethiopia.

The Tigray region, where we found around eight per cent of the indigenous cattle population of Ethiopia, is considered as the historic centre of the country, with the ancient pre-Aksumite and Aksumite civilisations in contact with the civilisations of the Fertile Crescent and the Indian subcontinent. Here, we used whole genome sequencing data to characterise the genomic diversity, relatedness, and admixture of five cattle populations (Abergelle, Arado, Begait, Erob, and Raya) indigenous to the Tigray region of Ethiopia. We detected 28 to 29 million SNPs and 2.7 to 2.9 million indels in each population, of which 7% of SNPs and 34% of indels were novel. Functional annotation of the variants showed around 0.01% SNPs and 0.22%-0.27% indels in coding regions. Enrichment analysis of genes overlapping missense private SNPs revealed 20 significant GO terms and KEGG pathways that were shared by or specific to breeds. They included important genes associated with morphology (SCN4A, TAS1R2 and KCNG4), milk yield (GABRG1), meat quality (MMRN2, VWC2), feed efficiency (PCDH8 and SLC26A3), immune response (LAMC1, PCDH18, CELSR1, TLR6 and ITGA5), heat resistance (NPFFR1 and HTR7) and genes belonging to the olfactory gene family, which may be related to adaptation to harsh environments. Tigray indigenous cattle are very diverse. Their genome-wide average nucleotide diversity ranged from 0.0035 to 0.0036. The number of heterozygous SNPs was about 0.6-0.7 times higher than homozygous ones. The within-breed average number of ROHs ranged from 777.82 to 1000.45, with the average sum of the length of ROHs ranging from 122.01 Mbp to 163.88 Mbp. The genomic inbreeding coefficients differed among animals and breeds, reaching up to 10% in some Begait and Raya animals. Tigray indigenous cattle shared a common ancestry with Asian indicine (85.6%-88.7%) and African taurine (11.3%-14.1%) cattle, with very small, if any, European taurine introgression. This study identified high within-breed genetic diversity representing an opportunity for breeding improvement programs and, also, significant novel variants that could increase the number of known cattle variants, an important contribution to the knowledge of domestic cattle genetic diversity.

A structural variant at ASIP associated with the darkness of hair coat is found in multiple zebu cattle populations.

A structural variant (SV) in the agouti signaling protein gene (ASIP), namely ASIP-SV1, has been found to strongly correlate with the darkness of hair coat in specific regions of the body of bulls from the zebu (Bos indicus) Nellore breed. Here we visually analyzed the whole-genome sequence of zebu and taurine (Bos taurus) cattle to elucidate the extent of spread of ASIP-SV1 in different cattle populations. Of 216 sequences analyzed, 63 zebu (45.9%) and five taurine (6.3%) animals had at least one copy of ASIP-SV1. Four of the taurine animals presenting the SV were Romagnola cattle, a breed with history of zebu introgression. The remaining taurine animal was a Simmental, a breed frequently used in crossbreeding. These data provide evidence that ASIP-SV1 is commonly found in zebu populations, in addition to taurine animals with zebu admixture.

Yolk Fatty Acid Content, Lipid Health Indices, and Oxidative Stability in Eggs of Slow-Growing Sasso Chickens Fed on Flaxseed Supplemented with Plant Polyphenol Extracts.

Previous attempts to increase the level of flaxseed in hens' diet for the production of n-3 polyunsaturated fatty acids (n-3 PUFAs)-enriched eggs have been commonly associated with undesirable effects on production efficiency, lipid health indices, and oxidative stability of eggs, requiring adequate research attention. This study investigated the effects of feeding a moderate level of flaxseed (FS) and plant polyphenol extracts (PPEs) on fatty acid content, oxidative stability, and lipid health indices in eggs of slow-growing Sasso T451A laying hens. One hundred and five hens were assigned to five groups (seven replicates of three) and fed on FS (75 g flaxseed and no antioxidants), VE8 (75 g flaxseed and 800 mg vitamin E), TS8 (75 g flaxseed and 800 mg Thymus schimperi), DA8 (75 g flaxseed and 800 mg Dodonaea angustifolia), and CD8 (75 g flaxseed and 800 mg Curcuma domestica) extract per kg diets. The egg yolk content of eicosapentaenoic acid (EPA, C20:5 n-3) in the DA8, TS8, and CD8 diets and docosahexaenoic acid (DHA, C22:6 n-3) in TS8 and CD8 diets significantly (p < 0.05) increased compared with the FS diet. The FS diet significantly increased the malondialdehyde (MDA) content in egg yolks, whereas the TS8 diet decreased it by 67% (p < 0.05). Little difference was observed in yolk fatty acid content between cooked and raw eggs. Production of n-3 PUFA-enriched eggs with favorable lipid health indices was possible through inclusion of PPEs extracted from local plant species grown in Ethiopia and a moderate dose of flaxseed in the diet of laying hens.

Genetic variation and demographic history of Sudan desert sheep reveal two diversified lineages.

More than 400 million sheep are raised on the African continent, the majority of which are indigenous and are primarily reared for sustenance. They have effectively adapted to various climatic and production environments, surviving and flourishing. The genetic relationships among these sheep populations remain understudied. Herein, we sequenced the entire mitochondrial DNA control region of 120 animals from Hamary and Kabashi and their crossbreed (Hamary x Kabashi) of Sudan desert sheep (SDS) to understand their maternal-inherited genetic variation and demographic history profiles and relate those to the history of sheep pastoralism on the African continent. The results show a diversified and predominant D- loop haplogroup B (n = 102, 85%), with all other sequences belonging to haplogroup A. Most of the maternal genetic variation was partitioned between haplogroup (76.3%) while within haplogroup accounted for 23.7% of the variation. However, little genetic differentiation was observed among the two breeds and their crosses, with our results supporting a Hamari maternal origin for the crossbreed. Bayesian coalescent-based analysis reveals distinct demographic history between the two haplogroups, two breeds and their crosses. Comparison of the two haplogroup showed that haplogroup B experienced an earlier expansion than haplogroup A. Unlike the breed-based comparison, the expansion of the two breeds started roughly at the same time, around 6500 years ago, with Kabashi having a slightly greater effective population size. The maternal ancestors of SDS may have diverged before their introduction to the African continent. This study provides novel insights into the early history of these two main breeds of Sudan desert sheep and their crosses.

The usefulness of maternally inherited genetic markers for phylogeographic studies in village chicken.

Phylogeography plays a major role in understanding micro and macroevolutionary processes dealing with evolutionary interpretations of geographical distribution. This field integrates information from molecular genetics, population genetics, demography, and phylogeny for the interpretation of the geographical distribution of lineages. The full mtDNA sequence and W chromosome polymorphisms were exploited to assess the usefulness of two maternally-inherited genetic markers for phylogeographic studies of village chickens. We studied 243 full mtDNA sequences from three countries (Iraq, n = 27; Ethiopia, n = 211; and Saudi Arabia, n = 5) and a 13-kb fragment of the W chromosome from 20 Iraqi and 137 Ethiopian female chickens. The results show a high level of genetic diversity for the mtDNA within and among countries as well as within populations. On the other hand, sequence analysis of the W chromosome shows low genetic diversity both within and among populations. Six full mtDNA haplogroups (A, B, C1, C2, D1, and E1) were observed and 25 distinct W haplotypes. The results support the effectiveness of full mtDNA sequences but not the W chromosome in tracing the maternal historical genome background with, however, weak within a country phylogeographic signal.

MicroRNA expression profile of chicken jejunum in different time points Eimeria maxima infection.

Coccidiosis stands as a protozoan disease of notable economic impact, characterized by an intracellular parasite that exerts substantial influence over poultry production. This invasion disrupts the integrity of the enteric mucosa, leading to the emergence of severe lesions and diminishing the efficiency of feed utilization in chickens. MicroRNA (miRNA) are short, non-coding RNA molecules with approximately 21-24 nucleotides long in size that play essential roles in various infectious diseases and inflammatory responses. However, the miRNA's expression patterns and roles in the context of Eimeria maxima infection of chicken intestines remain unclear. miRNA sequencing was employed to assess the miRNA expression profile in chicken jejunum during E. maxima infection. In this study, we analyzed miRNA expression profiles related to the host's immune response in the chicken jejunum during E. maxima infection. At 4 days infection and control (J4I versus J4C), 21 differentially expressed miRNAs in the jejunum were identified, comprising 9 upregulated and 12 downregulated miRNAs. Furthermore, in the jejunum, at 7 days infection and control (J7I versus J7C) groups, a total of 35 significantly differentially expressed miRNAs were observed, with 13 upregulated and 22 downregulated miRNAs. The regulatory networks were constructed between differentially expressed miRNA and mRNAs to offer insight into the interaction mechanisms between chickens and E. maxima coccidian infection. Furthermore, within the comparison group, we obtained 946, 897, and 281 GO terms that exhibited significant enrichment associated with host immunity in the following scenarios, J4I vs. J4C, J7I vs. J7C, and J4I vs. J7I, respectively. The KEGG pathway analysis indicated notable enrichment of differentially expressed miRNAs in the jejunum, particularly in J4I vs. J4C; enriched pathways included metabolic pathways, endocytosis, MAPK signaling pathway, regulation of actin cytoskeleton, and cytokine-cytokine receptor interaction. Moreover, in J7I vs. J7C, the KEGG pathway was significantly enriched, including metabolic pathways, protein processing in the endoplasmic reticulum, ubiquitin-mediated proteolysis, and FoxO signaling pathway. A comprehensive understanding of the host genetic basis of resistance with a combination of time-dependent infection to the Eimeria parasite is crucial for pinpointing resistance biomarkers for poultry production.

Genomic adaptation of Ethiopian indigenous cattle to high altitude.

The mountainous areas of Ethiopia represent one of the most extreme environmental challenges in Africa faced by humans and other inhabitants. Selection for high-altitude adaptation is expected to have imprinted the genomes of livestock living in these areas. Here we assess the genomic signatures of positive selection for high altitude adaptation in three cattle populations from the Ethiopian mountainous areas (Semien, Choke, and Bale mountains) compared to three Ethiopian lowland cattle populations (Afar, Ogaden, and Boran), using whole-genome resequencing and three genome scan approaches for signature of selection (iHS, XP-CLR, and PBS). We identified several candidate selection signature regions and several high-altitude adaptation genes. These include genes such as ITPR2, MB, and ARNT previously reported in the human population inhabiting the Ethiopian highlands. Furthermore, we present evidence of strong selection and high divergence between Ethiopian high- and low-altitude cattle populations at three new candidate genes (CLCA2, SLC26A2, and CBFA2T3), putatively linked to high-altitude adaptation in cattle. Our findings provide possible examples of convergent selection between cattle and humans as well as unique African cattle signature to the challenges of living in the Ethiopian mountainous regions.

Transcriptome analysis of differentially expressed circRNAs miRNAs and mRNAs during the challenge of coccidiosis.

Avian coccidiosis is a common enzootic disease caused by infection of Eimeria species parasites. It causes huge economic losses in the global poultry industry. Current control using anticoccidial drugs or vaccination is limited due to drug resistance and the relatively high cost of vaccines. Improving host genetic resistance to Eimeria species is considered an effective strategy for improved control of coccidiosis. Circular RNAs (circRNAs) have been found to function as biomarkers or diagnoses of various kinds of diseases. The molecular biological functions of circRNAs, miRNAs, and mRNAs related to Sasso chicken have not yet been described during Eimeria species challenge. In this study, RNA-seq was used to profile the expression pattern of circRNAs, miRNAs, and mRNAs in spleens from Eimeria tenella-infected and non-infected commercial dual-purpose Sasso T445 breed chickens. Results showed a total of 40 differentially expressed circRNAs (DEcircRNAs), 31 differentially expressed miRNAs (DEmiRNAs), and 820 differentially expressed genes (DEmRNAs) between infected and non-infected chickens. Regulatory networks were constructed between differentially expressed circRNAs, miRNAs, and mRNAs to offer insights into the interaction mechanisms between chickens and Eimeria spp. Functional validation of a significantly differentially expressed circRNA, circMGAT5, revealed that circMGAT5 could sponge miR-132c-5p to promote the expression of the miR-132c-5p target gene monocyte to macrophage differentiation-associated (MMD) during the infection of E. tenella sporozoites or LPS stimulation. Pathologically, knockdown of circMGAT5 significantly upregulated the expression of macrophage surface markers and the macrophage activation marker, F4/80 and MHC-II, which indicated that circMGAT5 might inhibit the activation of macrophage. miR-132c-5p markedly facilitated the expression of F4/80 and MHC-II while circMGAT5 could attenuate the increase of F4/80 and MHC-II induced by miR-132c-5p, indicating that circMGAT5 exhibited function through the circMGAT5-miR-132c-5p-MMD axis. Together, our results indicate that circRNAs exhibit their resistance or susceptive roles during E. tenella infection. Among these, circMGAT5 may inhibit the activation of macrophages through the circMGAT5-miR-132c-5p-MMD axis to participate in the immune response induced by Eimeria infection.

Ecological niche modelling as a tool to identify candidate indigenous chicken ecotypes of Tigray (Ethiopia).

The Tigray region is an ancient entry route for the domestic chickens into Africa. The oldest African chicken bones were found in this region at Mezber, a pre-Aksumite rural farming settlement. They were dated to around 800-400 BCE. Since then, the farming communities of the region have integrated chicken into their livelihoods. The region is also recognised for its high chicken-to-human population ratio and diverse and complex geography, ranging from 500 to 4,000 m above sea level (m.a.s.l.). More than 15 agro-ecological zones have been described. Following exotic chicken introductions, the proportion of indigenous chicken is now 70% only in the region. It calls for the characterisation of indigenous Tigrayan chicken ecotypes and their habitats. This study reports an Ecological Niche Modelling using MaxEnt to characterise the habitats of 16 indigenous village chicken populations of Tigray. A total of 34 ecological and landscape variables: climatic (22), soil (eight), vegetation, and land cover (four), were included. We applied Principal Component Analysis correlation, and MaxentVariableSelection procedures to select the most contributing and uncorrelated variables. The selected variables were three climatic (bio5 = maximum temperature of the warmest month, bio8 = mean temperature of the wettest quarter, bio13 = precipitation of the wettest month), three vegetation and land cover (grassland, forest land, and cultivated land proportional areas), and one soil (clay content). Following our analysis, we identified four main chicken agro-ecologies defining four candidates indigenous Tigrayan chicken ecotypes. The study provides baseline information for phenotypic and genetic characterisation as well as conservation interventions of indigenous Tigrayan chickens.

Influence of Eimeria maxima coccidia infection on gut microbiome diversity and composition of the jejunum and cecum of indigenous chicken.

Coccidiosis is an economically significant protozoan disease and an intracellular parasite that significantly impacts poultry production. The gastrointestinal tract microbiota plays a central role in host health and metabolism, and these microbes enhance chickens' immune systems and nutrient absorption. In this study, we analyzed the abundance and diversity of microbiota of the jejunum and cecum of a dual-purpose indigenous Horro chicken following Eimeria maxima infection. We compared microbial abundance, composition, and diversity at the 4- and 7- days post-infection using 16S rRNA gene sequencing. We obtained, on average, 147,742 and 132,986 high-quality sequences per sample for jejunum and cecum content, respectively. Firmicutes, Proteobacteria, Campilobacterota and Bacteroidota were the major microbial phylum detected in the jejunum content. Firmicutes were the dominant phylum for 4- and 7-days jejunum control groups accounting for (>60% of the sequences). In the infected group Campilobacterota was the dominant phylum in the jejunum (> 24% of sequences) at 4-and 7-days post-infection groups, while Proteobacteria was predominant at 4- and 7-days post-infection of the cecum (> 40% of the sequences). The microbial genus Lactobacillus and Helicobacter were found in the jejunum, while Alistipes, Barnesiella and Faecalibacterium were detected in the cecum. In the jejunum, Helicobacter was dominant at 4 -and-7 days post-infection (≥24%), and Lactobacillus was dominant at 4 -and 7- days in the control group (> 50%). In 4- and 7-days post-infection, Alistipes genus was the more prevalent (> 38%) in the cecum. Thus, clear differences were observed in the bacterial microbiota distribution and abundance between the jejunum and cecum, as well as between infected and control groups for both tissues. The results indicate that chicken intestinal microbial imbalance (dysbiosis) is associated with Eimeria parasite infection and will likely affect the host-microbial non-pathogenic and pathogenic molecular interactions.

Natural adaptation and human selection of northeast African sheep genomes.

African sheep manifest diverse but distinct physio-anatomical traits, which are the outcomes of natural- and human-driven selection. Here, we generated 34.8 million variants from 150 indigenous northeast African sheep genomes sequenced at an average depth of ∼54× for 130 samples (Ethiopia, Libya) and ∼20× for 20 samples (Sudan). These represented sheep from diverse environments, tail morphology and post-Neolithic introductions to Africa. Phylogenetic and model-based admixture analysis provided evidence of four genetic groups corresponding to altitudinal geographic origins, tail morphotypes and possible historical introduction and dispersal of the species into and across the continent. Running admixture at higher levels of K (6 ≤ K ≤ 25), revealed cryptic levels of genome intermixing as well as distinct genetic backgrounds in some populations. Comparative genomic analysis identified targets of selection that spanned conserved haplotype structures overlapping clusters of genes and gene families. These were related to hypoxia responses, ear morphology, caudal vertebrae and tail skeleton length, and tail fat-depot structures. Our findings provide novel insights underpinning morphological variation and response to human-driven selection and environmental adaptation in African indigenous sheep.